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1.  Introduction 


We  study  in  this  paper  the  propagation  of  an  electromagnetic 

wave  into  an  unbounded  domain  occupied  by  a  nonlinear  dielectric 

substance.  In  some  Cartesian  coordinate  system  (x^)  on  the  domain 
3 

S'  i  R  we  assume  that  the  domain  consists  of  the  half-space 
x^  >  0 ;  the  propagating  electromagnetic  wave  is  assumed  to  be  of 
the  form 

(1.1)  E  =  (0,E2(xl5t) ,0) ,  =  (0,0,H3(x1,t)  ) 

and  the  properties  of  the  nonlinear  dielectric  substance  occupying 
Si  are  delineated  by  the  nonlinear  constitutive  relations 


(1-2)  D  =  &(£)£,  B  =  pQpH 

which  determine,  respectively,  the  electric  displacement  and  magnetic 
field  in  S2  in  terms  of  the  electric  field  and  magnetic  intensity. 
The  constitutive  relations  (1.2)  are  commonly  considered  in  work  on 
nonlinear  optics  [1],  [2]  with  e  >  0,  p  >  0  being  scalar-valued 
functions  of  their  respective  fields.  As  £  =  e^E  +  £(£)  where 
£  is  the  macroscopic  polarization  vector  and  is  the  per¬ 

mittivity  of  free  space,  and  £(£)  =  X(£)£,  where  x(£)  is  the 
susceptibility,  £  =  (eQ+X(£))£.  For  an  isotropic  material, 

X  =  X(||E||),  p  =  |jl( ffjtjll )  and,  in  most  treatments  in  nonlinear  optics 
texts  these  functions  are  expanded  in  series  of  the  form 


(1.3) 


X0  +  xif:  + 


+  pH  +  0(H2),  H 
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where  Pg  is  the  permeability  of  free  space,  Xg  is  the  linear 

susceptibility,  X.^  the  first  nonlinear  susceptibility,  and  so  forth. 

\ 

Experimentally  X^  «  X^  so  that  the  presence  of  the  term  of  first 
order  in  E  in  the  expression  for  X  is  of  most  interest  when  E 
is  large  (e.g.  in  a  laser  beam)  while  for  most  purposes  one  may 
assume  that  p  =  p^ .  Many  of  our  results  apply  to  materials  in  which 
P (JtJ )  and  e(£)  =  6g  +  X(£)  are  more  general  than  the  forms  implied 
by  (1.3). 

In  addition  to  the  forms  (1.1),  (1.2),  respectively,  for  the 
electromagnetic  wave  propagating  in  S2,  and  the  constitutive  relations 
which  delineate  the  dielectric  material  that  occupies  ft,  we  need 
Maxwell's  equations  which  we  take  in  the  form 

-V  x  £,  div  £  =  0 

V  x  div  £  =  0 

We  also  introduce  the  notation 

e(£)  =  e((0,S,0)),  p(£)  =  |t((0,0,g)),  y,  €  R1 

and  assume  that 

h(i)  7(0,  7(0  €  C1(R1) 

h(ii)  7(0  t  0,  7<0  *  0,  V£  €  R1 

h(iii)  (s7(0)'  >  0.  (sp(0>’  >  0,  Vs  €  R1 

K  with 


or,  at  least,  for  all 
|^,|  sufficiently  small 


r  a* 

at 


(1.4) 


dD 

It 


3 


From  (1.2)  and  the  assumed  form  of  the  wave  (1.1)  we  obtain 


D  =  (0,D2(x1,t>,0),  B  =  (0,0,B3(x1,t)) 


with 


(1.5) 


D2 ( X! ’ t )  =  e(E^(x1,t))E2(x1,t) 


B3(x1,t)  =  |i(H3<x1,t))H3(x1,t) 


By  hypotheses  h(III)  these  relations  may  be  inverted  so  as  to  yield 


(1.6) 


E„(x  ,t)  =  - -  D  (x  ,t)  =  X(D  (x  ,t) )D  (x  ,t) 

2  1  e(E<D2(x1,t))  21  21  21 


H„(x,,t)  =  - -  B-(x, ,t))  =  y(Bq(x,  ,  t))B.,  (x.  ,t) 

3  1  n(H(B3(x1t))  31  31  31 


where  6  Rx  we  have  defined 


ms)  =  c -  ,  r(s)  =  z — = — 

e(E(£) )  p(H(^)) 


with  (for  p,  X,  €  R  ) 


P  =  e(£)£  -*•  x.  =  E(  p)  =  E(e(^)£)  =  X(p), 


p  =  •+  K  =  Mp)  =  H(p(^)x.)  £  Y (p)| 


Therefore  Vp  i  R1  (or  Vf  wi  th  |  #»  |  sufficiently  small) 


(1.7)  (pX<p))’  £  =  ~ 


d K  _  dE(P)  _  1 

dp  dp  ~  dp/dK, 


4 


by  h(iii)  and,  similarly,  (py(p))'  >0,  Vp  (  R1.  Note  that  in  the 
special  case  where  p  =  ,  y  =  1/Vg  >  0,  In  the  next  section  we 

demonstrate  that  shock  waves  may  be  expected  to  form  in  the  electro¬ 
magnetic  wave  (1.1)  as  it  propagates  into  the  half-space  >  0; 

We  also  derive  some  estimates  for  the  maximal  time  of  existence  of  a 
C1  wave  and  estimate  the  distance  travelled  by  the  wave  into  the 
half-space  until  the  development  of  the  shock;  these  latter  results 
are  presented  in  §3. 

There  is  considerable  literature,  ref.  [10]-[17],  on  shock  develop¬ 
ment  in  electromagnetic  theory  with  the  papers  that  are  closest  to  the 
present  work,  in  spirit,  being  those  of  Broer  [10],  [11],  Katayev  [16], 
and  Jeffrey  [11],  [12].  The  work  of  Broer,  however,  is  not  applicable 
to  those  important  situations  in  which  the  quasilinear  evolution  equa¬ 
tions  are  not  genuinely  nonlinear  while  our  work  differs  from  that  of 
Jeffrey  by  virtue  of  the  fact  that  by  working  with  the  constitutive 
relations  (1.2)  and  employing  D  and  B  as  our  basic  variables, 
instead  of  E  and  B,  as  in  [12],  or  E  and  H,  as  in  [11],  we  are 
able  to  write  our  evolution  system  in  conservation  form;  the  importance 
of  working  with  D  instead  of  E  was  emphasized  in  [6]  and  is  based 
on  the  fact  that  E  is ,  in  general,  not  divergence  free  in  a  nonlinear 
dielectric  while  D  is  if  there  is  zero  free  charge.  Some  of  the  con¬ 
sequences  of  having  the  equations  for  D,  B  in  conservation  form, 
i.e.,  the  implications  of  the  Rankine-Hugoniot  and  Lax  k-shock  condition 
are  developed  in  the  nc?xt  section.  The  application  of  Lax’s  elegant 
work  [3],  also  simplifies  the  asymptotic  estimate  for*  t  in  <j?  a?  I 

makes  possible  a  rather  simple  and  explicit  computation  in  §3  S 

the  distance  travelled  by  the  wave  into  the  half-space  before  shock 


deveJ opment  occurs . 


■rtiii 


2 .  Shock  Development  and  Propagation 

To  simplify  the  notation  in  this  section  we  set  x  =  x^, 

D  =  ,  E  =  H  =  H.^,  B  =  In  view  of  the  forms  of  the 

electromagnetic  field  vectors  in  the  wave  entering  the  half  space 
x^  >  0  Maxwell's  equations  (1.4)  reduce  to  the  pair  of  equations 


f9  d_D  3_E  _  aH  .  3_B  aH 

aE  *  at  ‘  'ax  ’  aH  *  at 


3E 

3x 


or 


(2.2) 


(e(E)E)'Et  +  Hx  =  0 


(n(H)H)'Ht  +  Ex  =  0 


Setting,  a(^)  = 


>  0,  b(£)  = 


0 ,  €  R  , 


(e(£)^)'  (n(^)O' 

we  see  that  E(x,t),  H(x,t)  satisfy  the  first-order  quasilinear 


system 


(2.3) 


E .  +  a(E)H 
t  x 


H .  +  b(H)E  =  0 

t  x 


which  is,  unfortunately,  not  in  the  usual  conservation  form.  We 
thus  rewrite  the  system  in  the  form 


.TTeT  +  Hx 


0, 


Ht 

binr  +  r*x 


and  note  that  in  view  of  the  definitions  of  a(*),  b( • ) ,  and  (l.b) 


With  (2.7)  we  associate  initial  data  of  the  form 


1 
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(2.8)  D( x , 0 )  =  Dq(x),  B(x ,0 )  =  Bq(x) 

and  real  characteristics  in  the  x,t  plane 


(29)  —  -  +  E 1 ( D ( x , t ) ) 

dt  "  "  •/  ^0 

Note  that  for  £  e  .  E'(£)  =  (£\(£))’  >0  by  (1.7).  The  positivity 

of  E'(£),  V£  €  Rd  is  equivalent  to  the  strict  hyperbolicity  of  the 

system  (2.7);  via  standard  a  priori  estimates  (Lax  [3],  Nishida  [4])  on 

the  Riemann  Invariants  associated  with  the  system  (2.7)  (these  being 

defined  below)  it  can  be  shown  that  if  sup|Dn(x)|  and  sup(Bn(x)| 

1  u  1  u 

R  R 

are  sufficiently  small,  and  E'(0)  >  0,  then  for  as  long  as  a 

sufficiently  smooth  solution  (D(x,t),  B(x,t))  of  (2.7),  (2.8)  exists, 

on  0  <  t  <  t  for  instance,  we  will  have  E'(D(x,t))  >  0, 

max 

x  €  r\  0  1  t  <  t  .  For  sufficiently  small  initial  data,  therefore, 
the  (real)  characteristics  (2.9)  are  well  defined  on  the  maximal 
interval  of  existence  of  Cd  solution  even  if  only  local  hyperboli¬ 
city  of  (2.7)  obtains,  i.e.,  even  if  we  only  have  E'(0)  >  0. 

We  now  define  the  Riemann  Invariants  associated  with  (2.7)  to  be 


f  r(D,B) 


(2.10)  I 

s  ( D ,  B ) 


rD 

0 


VF’T^T  d£ 


si  ,ln 


v/PTST  dK 


By  standard  results  r  and  s  are  constant  along  their  respective 
characteristics ,  i.e., 


r'(x,t)  ,  S  ££  = 


ax 


(2.11)  ( 


s'  ( x , t )  =  —  +  —  (D(x,t~TT  as  _ 


at 


7 


lo 


ax 


and  we  have  the  following  blow-up  results  which  are  consequences, 
respectively,  of  the  work  of  Lax  [3]  and  Klainerman  and  Majda  l5J. 

(A)  If  Dq(x)  is  periodic  on  R1 ,  BQ(x)  5  and  £"(0)  *  0 

(so  that  the  problem  exhibits  genuine  nonlinearity)  then  finite-time 
blow-up  must  occur  for 


r(x,t)  =  B  ( x , t )  +  —  VE * (D(x , t ) )  D  ( x , t ) 

x  x  t —  x 

#0 


’•Vt 


D  (x,t)  +  —  7E  * ( D(x , t ) )  D  (x,t) 


-*■  7,  .  >D  s  (D  ,D  )  must  blow-up  in  finite  time  and  a  shock  develops 

^  X  ,  t  /  L  X 

Futhermore,  it  is  a  consequence  of  the  work  in  [3]  that 


(2.12)  t 


VE'(O) 


max  max  |  Dq  Cx7  j  |  E"  (  0 )  | 


a  result  we  will  apply  in  §3  to  a  specific  class  of  nonlinear  dielec¬ 
tric  materials. 

(ft)  rjuppor."  that  l)f (x),  B  (x)  both  have  compact  support  in  R*  ; 
then  so  will 


I 
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r(x ,0 ) 


s (x  ,  0 ) 


B„(x)  *  —  f  0  i\K 

J(> 


B  (x)  -  —  fDo<I!>  ,/E'ir. )  dr. 

A  J° 


From  the  recent  work  of  Klainerman  and  Kajda  [5]  it  then  follows  that 
if  r(#  ,0) ,  s(*,0)  are  also  of  class  then  any  solution  of 

the  initial-value  problem  for  the  diagonalized  system  (2.11)  must 
develop  singularities  in  finite  time  in  the  first  derivatives 
r  ,  if  Ef(£)  is  not  constant  on  any  open  interval.  For  example, 

if  with  >  0,  X(£)  =  Xq  +  X^£^  then  E(£)  =  *  Clearly 

E  f  (  0  )  =  Xq  >  0  (local  hyperbolic!  ty )  and  E 11  (  0  )  =  0  (loss  of  genuine 
nonlinearity)  but  EM(£)  =  if  C  /  0,  so  that  the  result 

of  [5]  applies  for  ,  compactly  supported  initial  data. 

Remarks .  If  (D,B)  is  a  sufficiently  smooth  solution  of  (2.7), 

2 

say  class  C  ,  then  clearly  we  may  eliminate  so  as  to  obtain  the 
scalar  nonlinear  wave  equation  for  D(x,t): 


(2.13)  D(  x ,  t ) 

dt 


1/H  ■  E(D(x,t)) 

dX 


This  equation  was  derived  by  this  author  in  [6]  by  specializing  the 
three  dimensional  evolution  equations 


(2.19) 


o?r>. 

l° 


V  (  X (£>) ) -grad ^( grad X(D)  *D)  , 


obtained  under  the  assumption  that  £  =  £  =  in  &  •> 


to 
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the  casio  of  ari  electromagnetic  wave  ot  the  form  (1.1)  propagating  through 
an  nonlinear*  dielectric  cylinder,  with  the  direction  of  propagation 

directed  along  the  axis  of  the  cylinder.  Several  facts  may  be  noted 
about  the  simple  wave  equation  (2.13) 

1)  Suppose  we  set  D(x,t)  =  Gx<x,t)  in  (2.13).  Then  provided  that 
G(x,t)  is  sufficiently  smooth  we  obtain 

Sttx  =  l,»0«Gx>xx  " 

Gtt  *  l/“oE«Vx  *  K(t) 


whet'e  K(t)  is  an  arbitrary  function  of  t.  Thus  the  usual  scalar 
nonlinear  wave  equation  which  arises  in  one-dimens ional  motions  of  a 
nonlinear  elastic  body  is  not  equivalent  to  (2.13). 

2)  Suppose  that  D(x,t)  is  a  solution  of  class  of  an  initial- 

value  problem  on  associated  with  (2.13)  and  is  such  that  for 


(2.15) 


E ( D( -« 


t ) )  dt 
x 


< 


00  where 


E(D(  —  ,t)) 

x 


lint  E( D(x,t))] 

X-^-oo 


If  we  set 


(V.16)  §(x,t) 


"hr 


X 

Dt(y  ,t)dy 


»t 

F(D(  —  ,t:))  dt  +  C 
.  x 

*  —  ao 


where  C  is  an  arbitrary  constant,  then 


1] 


b  (  x  ,  t  )  =  -|j  ..  U .  (  x  ,  t )  -and  by  (  2 . 1  3  ) 

x  0  t  * 


B  (x,t ) 


x  .2 


f* 

-m-q  — 2  D(y»t)  dv  "  E(D(-*,t))x 

' at 


rx  2 

E(D*y,t))dy  -  E(D(  -<'• ,  t )  ), 

J  a*/ 


Qy 


E  (D(x  ,t )  ), 


so  that  the  pair’  (D(x,t),  B(x,t))  is  a  solution  of  the  system  (.2.1) 
This  equivalence  between  the  system  (2.7)  and  the  scalar  nonlinear 
wave  equation  (2.13)  only  holds  if  B  is  well-defined,  i.e.,  if 
(2.15)  obtains.  As  . 


E(D(-»,t))  =  lim  [E’(D(x,t))D  (x,t)], 

x  x^-»  x 


if  (2.13)  is  strictly  hyperbolic,  i.e.,  Ef(£)  >0,  €  R  ,  then 

clearly  (2.15)  will  obtain  if  lim  |D  (x,t)j  =0,  0  .  t  <  t 

J  .  1  x  ’  1  max 

We  now  return  to  the  system  (2.7);  shock  development  and  propagation 
in  the  more  general  system  (2.5),  which  is  strictly  hyperbolic 
provided  E'  (<)//'  (£)  >  0,  ^  ,  will  be  dealt  with  in  a  forth¬ 

coming  paper  [7],  If  we  let  s  denote  the  speed  of  the  shock  which 
develops  either  in  case  (A)  or  case  (B)  above,  and  [F]  the  jump  in  quan¬ 
tity  F  across  the  shock  then  the  Rankine-Hugoniot  conditions  require 
that  sfu^]  =  f fR  1 ,  k  =  1,2  where 


U  =  lD\ 
\B  I 


l/p0  B 
\(D)D 
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We  therefore  obtain  the  conditions 


(2.17)  < 

s  [  D] 

=  —  [B]  = 

0 

[H] 

s  [  B  ] 

V 

=  [X(D)D]  = 

[E] 

from  which  follow  the  simple 

relati 

ons 

(i) 

[D][E] 

=  [  H  ]  [  B  ]  = 

1 

^0 

[B]2 

( ii ) 

s2[D] 

It 

( — 1 

CD 

i — i 

mU° 

II 

1 

*0 

[E] 

The  shock  speed 

is  ,  therefore 

,  given  by 

(2.18) 

s  = 

♦  1  /CE] 

f«r  ^ [D] 

5  t 

1 

TOdTdT 

[D1 


so  that  two  shocks  are  possible,  one  moving  to  the  left  and  one 
moving  to  the  right.  We  now  apply  Lax’s  [8],  [9]  k-shock  conditions 

=  0,  i  =  1,2;  i.e.,  we  require  that 

for  either  k  =  1  or  k  =  2 


to  the  system  u.  +  —  f. 

St  i  ax  i 


(2.19)  \K(u_) 


where 


*  s  >  ( u+ ) 

u  =(£),  ?(u)  =  (1/|J-0B\ 

\x(d)d|. 


and  the  X,  ,  k  =  1,2  are  the 
k 


distinct  real  eigenvalues  of  the  matrix  However,  by  (2.9) 


(2.20) 


^  t/E  *  (  D)  , 
/“O 


/d  \ 

-► 

/  D.\ 

> 

\BJ’ 

u+ 

= 

(»:) 

-  ^  ^ftdT 

denote,  respectively,  the 
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values  of  u  behind  and  in  front  of  the  shock.  The  conditions 
(2.19)  represent  one  formulation  of  a  classical  entropy  condition 
for  solutions  of  hyperbolic  conservation  laws  containing  a  shock. 

Using  (2.20),  (2.19)  becomes 


(2.21), 


VE'iuj  >  s  >  \/etTd  ) 

-  VF* TFT  >  S  >  -  VFTd^T 


or 


However,  by  the  definition  of  E,  E  =  E(D)  =  E(D(E)),  where 
0(E)  =  f.(E)E  so  that  (2.21)  is  equivalent  co 


(2.22)/ 


/0  '  ( E  ) 


V0'(E  ) 


>  v^O  s  > 


\JV  'Te+~7 


or 


>  s  >  - 


V0 ' (E+  ) 


For  the  shock  moving  to  the  right,  sp  ~  ~Z  J  [  jflEyj 


so  that  only  the  first  inequality  in  (2.22)  makes  sense,  and  we 
must  have 

(2.23a)  V0 '  ( E_ )  <  <  V0’  (E+  ) 

For  the  shock  moving  to  the  left  s^  =  -  —  J~j 

so  that  only  the  second  inequality  in  (2.22)  makes  sense,  and  wo 
obtain 

(2.23b)  v/D tTtT)  >  /"VeT"  >  ^ '  ( E+  ) 


J 
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Thus  L  he  k -shook  conditions  oi  Lax  predict  that  (2.23a)  must  hold 
for  the  shock  moving  to  the  right  while  (2.23b)  must  hold  for  the 
shock  moving  to  the  left.  We  now  examine  the  implications  of  (2.23a,b) 
for  the  simple  but  physically  important  case  where 

(2.24)  e(E)  =  e o  +  E^  ,  >  0,  c ^  >  0 

We  recall  that  by  the  definitions  of  E,  V,  €  R  ,  %  =  E(P(£)) 

where  P(£)  =  7 (.OK  s  eQK  +  Thus  ViO  =  0  if  and  only  if 

K.  -  0.  A  direct  computation  yields 


so  that  E" ( P( 0  ) )  =  E" ( 0 )  =  0  but  E"(P(£))  *  0  for  all  K  *  0 . 
Thus  for  initial  data  with  compact  support  the  results  of 

Klainerman  and  Majda  [5],  as  previously  stated,  apply  and  shocks 
will  develop  in  finite  time.  As  P'(E)  =  +  3e^,E",  the  condition 

(2.32a)  relative  to  the  shock  moving  to  the  right  with  speed 


*v<vve;  +(e++j:_)e_>> 


(2.21)  ./tQ  +  3t2E2_  <  /e0  +  e2(Ej+(E++E_)E_) 


<  V  60+3e2E+ 


This  last  inequality  clearly  implies  that 


[  E2  +  E  E  > 


(2.28a) 


2  2 
E_  +  E+E_  <  2E^ 


from  which  we  easily  deduce  that 

/  n  o  n  t  \  o  r-i 2  t-i  2  r> 


(2.28b) 


E2  >  E  E  >  2E?  -  E2 

+  -  + 


2  2 

(2.28c)  E+  >  E  (across  the  shock  moving  to  the  right) 

In  a  completely  analogous  fashion  we  obtain,  for  the  shock  moving 
to  the  left 


"  t  =  -|Aq  ( eQ  +  e  2  (E^  +(E++E_)E_))  * 
while  across  this  shock 

( 2E l  >  E l  +  E+f_ 

(2.29a)  / 


2 

2  , 

2E_ 

>  E+  +  E+ 

2 

2E. 

<  r.'  +  e. 

from  which  we  deduce  that 


(2.29b)  2E7_  -  E7  >  E+E_  >  2E7  -  E^ 

and 

2  2 

(2.29c)  E_  >  E+  (across  the  shock  moving  to  the  left.) 

The  situation  corresponding  to  (2.28c)  and  (2.29c)  is  depicted  below 
where  we  have  denoted  the  shock  moving  to  the  right  with  speed 
by  x  =  x  (t)  and  the  shock  moving  to  the  left  with  speed  . 

L  v 

by  x  =  x  (t) 


_In  the  figure  sketched  above  E+ ,  E_  denote,  respectively,  the  values 
of  E(x,t)  in  front  of  and  behind  the  respective  shocks.  For  the 
shocks  xA(t),  *r(t)  therefore,  (2.28c)  and  (2.29c)  predict  that 
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'/ 

L‘  (x,0  must.  increase  as  we  cross  the  respective  shocks ,  moving  in 
the  direction  of  increasing  x.  However,  that  part  of  the  energy 
residing  in  the  electromagnetic  wave  which  depends  on  £  is  given 
by  = 

=  hV(L)E  =  ^((e0+e2E2)E)E 

and  thus  we  have  the  result  that 


(2.30)  ®r(x,t) 

k, 


^egE2(x,t)  +  (x  ,t ) 


must  increase  as  we  move  across  the  respective  shocks  in  the  direc¬ 
tion  of  increasing  x. 


Remarks.  For  the  equations  (2.3),  i.e. 


j  Et  +  a(E)Hx  =  0 

H  +  b(H)E  =  0 

t  x 


where  a(E)  >  0,  b(H)  >0  it  is  a  simple  matter  to  show  that  there 
exist  solutions  of  the  form  E(x,t)  =  ESx-Xt),  H(x,t)  =  H*(x-Xt) 
where  EQ(x)  =  E(x,0),  H0(x)  =  H(x,0).  In  fact  X  =  ±\/a(E)b(H) 

so  that  (implicit)  travelling  wave  solutions  of  the  form 


IE  =  Eg (x  ±  \/a(E)b(H)  t) 
H  =  HQ(x  ±  Va(E)b(H)  t) 


may  be  well-defined,  at  least  for  small  values  of  t .  If  the 

~  -1 

material  is  such  that  m- ( H )  =  |ig  then  b(H)  =  fig  and  (2.31) 

2  " 

ftolutjons  of  I  ho  form  (2.31),  for  l hr  system  (7.3)  h.ivr  discus 

by  Kataycv  L10|. 
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assumes  the  form 


(2.32) 


E  =  En ( x  ±  ^  ViTET  t) 


H  =  Hq(x  t  Up*5  \/a(E)  t) 


If  we  define 


F ( x , t , E )  =  E  -  En(x  ±  p"'5  VaTET  t) 


then  F(xn,0,0)  =  0,  for  any  value  of  xn  such  that  EQ(x0>  =  0, 


FE(x,t,E) 


1  ±  ^jXq^Eq  (xi^Q^VaTETt )  *(a(E))-!sa'(E) 


a(E) 


e(E)  +  c ' (E)E  -  a(0 ) 


(2.33) 


a '  ( F )  =  2e»(E)  +  e"(E)E  -  a'(0) 


2e' (0) . 


In  the  typical  situation  where  e(E)  =  eQ  +  e^E  ,  eQ  >  0,  c  9  >  0 , 
a(0)  =  f. 0  and  a'(0)  =  0  so  that  Fe(x0,0,0)  =  1.  If  e(E) 

=  e0  +  e1E  then  a(0)  =  eQ,  a'(0)  =  2r  and  Fe(x0,0,0) 

-V 

e  1^0 

=  1  ±  —  ■  Eq(Xq).  By  the  implicit  function  theorem,  therefore,  if 

Xg  is  such  that  Eq(x^)  =  0  and,  either. 


e0  +  e2E 


e0  +  elE’  and’  eo  4  ei(io  Eo'(xo)  *  u 


a  solution  E  =  E(x,t)  of  F(x,t,E)  =  0  will  exist  for 
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and  |t|  sufficiently  small. 

Suppose,  now,  that  we  differentiate  the  first  relation  in  (2.32) 
through  with  respect  to  x  and  solve  for  we  easily  obtain 


-1 


(2.34)  Ex  =  E^(x±p~VilE)t)/  l/iy  a' (E)t 


where  a(E),  a '(E)  are  given  by  (2.33).  If  3c*  >  0  such  that 


(2.35)  €  R1,  |^|  sufficient 

then  for  the  wave  moving  to  the  left  with  velocity 


iy 

M. 


0 


small 


';Va(E)' 


(2.36)  1  - 

as  t  -►  t*  =  2/y.o1  € '*  thus,  by  (2.34) 


(2.37)  E(x,t)*+  +  °°  as  t  -►  t* 

x 


It  therefore  follows  directly  from  (2.34)  that  a  shock  can  be  expected 
to  develop  provided  the  constitutive  relation  in  the  material  is  such 
that  (2.35)  is  satisfied.  For  the  case  in  which  e(£)  = 

£  €  r\  (2.35)  assumes  the  form  +  <  2e^/e*  which  is 

certainly  satisfied  for  |^|  sufficiently  small  if  c*  is  chosen 
sufficiently  small;  the  prediction  of  shock  development  via  (2.34) 
is  not  suprising  in  this  case  since  we  are  in  the  genuinely  nonlinear 
s i tua Lion . 

'v  2 

However,  if  c  (£)  -  then  (2.35)  is  equivalent  to 


/E0+3e2* 
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w'ni'.h  can  not  be  satisfied  7^.,  |^l|  sufficiently  small,  no  matter  how 

small  e*  >  0  is  chosen;  this  is,  of  course,  a  situation  in  which 

genuine  nonlinearity  fails  and  shock  development  can  not  be  shown  to 

follow  directly  from  (2.34).  This  is,  essentially,  the  sort  of 

situation  in  which  the  utility  of  the  result  proven  by  Klainerman 

and  Majda  [5]  becomes  apparent  provided  the  initial-value  problem, 

formulated  in  terms  of  the  fields  J3  and  J),  instead  of  E  and 

is  such  that  the  initial  data  are  and  compactly  supported. 

Analogous  results  may  be  presented  for  the  more  general  situation 

described  by  the  system  (2.3),  where  we  do  not  assume,  a  priori,  that 

2. 

p(£)  =  (  R  ,  but  a  discussion  of  such  results  will  be 

relegated  to  a  forthcoming  work  [7].  It  is  noteworthy  to  remark  that 
solutions  of  the  basic  form  (2.32),  for  an  electromagnetic  wave 
propagating  into  a  half-space  filled  by  a  nonlinear  dielectric  sub¬ 
stance  described  by  an  arbitrary  functional  relationship  between  D 
and  E,  were  previously  obtained  by  Broer  [10],  (17]  where  1  lie  analysis 
was  taken  as  far  as  the  determination  of  Ex  which,  it  was  shown, 
could  blow-up  in  finite-time  thus  leading  to  the  initiation  of  non¬ 
uniqueness  in  the  solution. 

Remarks .  If  we  consider  a  material  for  which  p(H)  =  pQ  then  (2.3) 
reduces  to 


which  is  easily  seen  to  imply  that  E(x,t)  satisfies  t  lie  nonlinear 
wave  equation 


2 2 


Some  Computational  Aspects  of  Shock  Development 

In  this  section  we  derive  some  estimates  for  t  ,  the  maximal 

max 

time  until  the  development  of  a  shock  in  an  electromagnetic  wave  of 
the  form  (1.1)  which  is  propagating  into  a  half-space  filled  with  a 
nonlinear  dielectric  substance.  We  also  estimate  the  velocity  of 
the  wave  in  the  nonlinear  dielectric  and  the  distance  travelled  by 
the  wave  into  the  half-space  until  the  shock  develops.  We  assume 
that  the  dielectric  conforms  to  the  constitutive  hypothesis 


D  “  e0E  +  P<E),  B  =  M-qH 

2 

with  P(E)  =  XgE  +  X^E  ,  Xq  >  0,  X^  >  0.  The  quantities 

X ^  have  been  defined  in  §1  as  being,  respectively,  the  linear'  and 

(first)  nonlinear  susceptibilities  of  the  dielectric.  Thus, 

D  =  r,  (E)E,  e(E)  =  +  e^E  with  eo~eO+X0’  c  t>  =  X  . 

2 

From  D  =  (^q+X^)E  +  X^E  we  easily  compute  that 


(3.1) 


E(x , t ) 


(VV 

2X , 


2X. 


J  (e0+X0)2+4X1D(x,t) 


Choosing  the  positive  sign  on  the  radical  and  expanding  in  a  power 

3 

series  we  obtain  E  =  \(D)D  +  0(D  )  where 

r 


(3.2) 


\(D) 


X0  ♦  X  D 


<Ko+xoi 


0,  \n 


VV3/2 


Our  problem,  up  to  terms  of  order  Lf, 


is  then  of  the  form  (2.7)  with 
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E<D)  =  "  '  ^W3'2  D? 


;o  that 


(3.3) 


X  i  2  X 

E'(D)  =  7—  t-y — <r  -  -7 - 77- r-3/o  D  >  0  for 

<VV  (eo+xo) 3/2 


I D  |  <  (e0+X0)"2/2X 


We  assume  that  initial  fields  of  the  form  (2.8)  are  prescribed  with 
Bq(x)  =  0,  Dq(x)  periodic  on  r\  and  sufficiently  small,  so  that 
the  standard  a  priori  estimates,  based  on  the  Riemann  invariants 
associated  with  (2.7),  imply  that  (3.3)  is  satisfied  for  as  lonp,  as 
a  C1  field  D(x,t)  exists.  Actually 


D0<*> 


<e0*X0>Et|(x)  *  Xg^x) 


so  that 


(3.4)  D0'(x) 


(WVx)  +  2XlE0(x)E0(x) 

e’(x)[(60+x0)+2x1E0(x)] 

2x1Eq(x)Eq(x) 


In  making  the  approximation  in  (3.4)  we  are  assuming  that  the 
initial  field  is  that  of  a  high  intensity  (laser)  beam  whose  strength 

g 

is  of  the  order  of  magnitude  10  (volts/meter)  while  the  per- 

—  1  3  ? 

jiiittivity  of  free  space  is  of  the  order  10  (Columbs/Newton-Met .err.  ) 

is  of  the  order  t 0~ 1  !  ( Co  I umbs /Newton-Meters ‘  ) ,  and  X(  is  of 

i  j  •>  -1  3 

the  order  10  (Columbs/vol  ts  ‘  ) .  For  example,  =  4  x  10  lor' 
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::  -.x  matched  KDP  [2].  The  various  units  (MKS  system)  are  related 

by  the  identification  Newtons/Columb  s  Volts/Meter,  which  are  the 
dimensions  of  the  electric  field  E  (that  of  D,  the  electric 
induction  field,  as  well  as  of  P,  the  polarization,  are  then 
Columbs/Meters2 ).  Thus,  in  (3.4),  the  quantity  U  +XQ)  is  of  order 
10~13  while  *2Eo  i s  order  10  4. 

Since  E"(0)  =  -2X 2 / ( c0+XQ ) 3/ 2  t  0  the  problem  (2.7),  (2.8) 
is  genuinely  nonlinear  and  the  results  of  §2  apply.  In  particular, 
we  find  that 


(3.5) 


t 

max 


M'0(f'0  +  X0) 


(max | EqEq | ) 


Also,  if  we  denote  by  v^  the  velocity  of  the  beam  in  the  dielectric 
then 


There  fore , 
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v/2^oX  1E 


^0X1 


(max  E) 


Finally,  the  distance  travelled  by  the  beam  into  the  dielectric 
ha: i -apace  until  the  shock  develops  is  given  approximately  by 


max 


max 


v'2m’0X1 


(max  E) 


-("5/2-0)  ^IVob"1 


/ 


a 


or 


v 


(5.7) 


max 


_1_ 


|io(eo+xo) 


2^ 

(max  E)  ^(maxlE^E^)) 


X 


-1 


\  b 

■  V 


,  -  'V 
\ 

'v7  >' 


\ 


It  over  the  distance  travelled  by  the  beam  into  the  ha# f -space, unt 1 1 
shock  development,  max  E^  max  E If  then  we  have  the  (admittedly  crude 

n  \ 

. .  ^  J 


buc,  nonetheless,  informative)  estimate 


( 


'  (3.8)  S  .. 

max  0 

,  ^  ^  f 


i — i  r  /  _  'j  /  /  “  a  _i 

Climax  E„)  (maxlE^j) 


V'1' 


0 


where 


0vo0  0 


is  a  characteristic  material  coefficient 


£/sy*l\  . U  i/v-C  2  f 


"  jv\/2  )  x: 

which  may  be  associated  with  a  particular  nonlinear  dielectric 

substance.  For  most  common  nonlinear  dielectric  substance,  e.g., 

index  matched  KDP-,  will  be  a  very  large  number,  something  of  the 

order  of  magnitude  10^},.  Thus  (1.8)  indicates  that  even  for  an  inc  i- 

9 

dent  high  intensity  beam  of  the  order  of  magnitude  ]0  vol  i  s/met er  a 
steep  gradient  on  the  incident  beam  will,  be  needed  so^lhat  shook 


development  may  occur  within  distances  .obtainable  under\laborvt  ory  ;  i 

(_  /if  U  r-<^  "f  \  .  it  '■ 

editions.  iC-  \  .  \  p  i  6 

t  .  .  V !  ,  i  ^  ■'tL  2.  f  r*  ,  ( 


£1^ 
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